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Abstract 

We provide, in an extremely simple way, an upper bound to the 
minimum number of unitary operators describing a general random- 
unitary channel. 



1 Introduction 

A channel £ — i. e. a completely positive trace-preserving map — acting on 
density matrices p defined on a finite dimensional input Hilbert space ^ 
(for sake of simplicity, we consider here channels with equal input and output 
Hilbert spaces; the generalization is straightforward) is called random-unitary 
if it admits a Kraus representation [T] as 

£{p) = Y.P^U,pUl (1) 

i 

where Pi are probabilities and Ui are unitary operators. This definition in- 
volves an existential quantifier, and there is no known constructive algorithm 
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to check whether a given channel is random- unitary or not. Only the nec- 
essary condition of being unital, that is, of preserving the identity matrix, 
^(1) = 1, holds!. 

Nonetheless, random-unitary channels play a very special physical role 
among all possible evolutions that an open quantum system can undergo [3] . 
In fact, Gregoratti and Werner |^ proved that they are the only irreversible 
channels that can be perfectly corrected using, as the only side- resource, 
classical information extracted from the environment. This property is ac- 
tually sufficient and necessary for a channel to be random-unitary and one 
would prefer to adopt this one as the physical and operational definition of 
random-unitary channels. The idea of using the environment as a resource 
then initiated investigations about environment- assisted capacities for quan- 
tum channels 00 13- R^f. [F the problem also of quantifying the amount 
of classical information needed to perfectly correct a random-unitary chan- 
nel was raised for the first time in the case of decohering evolutions. In 
fact, for a given random- unitary channel, the form is highly non-unique 
and the Shannon entropy H{pi) of the probability distribution weighing the 
unitaries Ui can be "artificially" made as large as desired. Consequently, in 
order to derive sensible information-theoretic relations regarding the infor- 
mation dynamics in a random-unitary evolution, one has to single out the 
random-unitary Kraus representation minimizing H{pi). 

In the present paper we derive an upper bound for the minimal num- 
ber of unitary operators needed in Eq. (^Q), thus providing also a bound to 
the amount H{pi) of classical information needed to be extracted from the 
environment in order to invert the random-unitary evolution. Our bound, 
proved for generic dimension, does not catch the peculiar geometry that bis- 
tochastic qubit channels enjoy: For d = 2, it is provably non tight. However, 
the qubit case is completely understood and all random- unitary qubit maps 
have already been explicitly characterized (see, e. g. Ref. 0). In this sense, 
a bound for the qubit case is completely useless. On the contrary, as soon 
as one leaves the two-dimensional world, already for d = 3, the bound we 
provide is generally non trivial. 

^For two-dimensional systems, a channel is random- unitary if and only if it is unital. 
For higher dimensional systems, if a channel is random- unitary it is also unital, but the 
converse does not hold |2j. 
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2 Properties of random-unitary channels 



Let us given a channel 8 acting on density matrices p defined on the input 
Hilbert space J^. As a consequence of the Stinespring theorem jTU], we can 
write it as follows [TT] 

£(p) = Tr4t/(p®|0)(0U)f/t], (2) 

namely, as a unitary interaction between the system and an ancilla (or en- 
vironment, described by the Hilbert space J^a)i followed by a trace over the 
ancillary degrees of freedom. If the ancilla input state is a pure one — like 
in Eq. Q — Gregoratti and Werner Q proved that, for all possible unitary 
interactions U in Eq. (j21), and for all possible decompositions of the channel 
£ into pure Kraus representations £{p) = ^^EipE], there exists a suitable 
rank-one POVM on the ancilla, let us call it {|'yj)(fi|a}, Yli kj)('^i|a = la, 
such that 

E;pE\ = Tr4^(p® |0)(OU)t/t (1 ® \v,){v,\a)]. (3) 

As an immediate consequence, if the channel £ admits a random- unitary 
decomposition as £{p) = ^iPiUipUj , with Ui unitary operators, there ex- 
ists a rank-one POVM on the ancilla, {|aj)(aj|a}, such that the probability 
distribution of its outcomes does not depend on the input state p, since 

Tr[t/(p® |0)(OU)f/t (l®|a,)(a4)] = Tr[p,t/,pt//]=p„ Vp,V2. (4) 

It is now useful to introduce the channel £ from density matrices on 
to density matrices on defined as 

^(p) = Tr^[f/(p®|0)(OU)t/t]. (5) 

Since the unitary interaction U is unique up to local isometries on J^, we can 
consider such an ancillary (or complementary channel as a canonical 
one. In turn, the channel £ acting on density matrices, induces a unique dual 
ancillary channel £' acting on operators Oa on M'a as follows 

Tr[m)p]=Tr[0„^(p)]. (6) 

This is nothing but the Heisenberg picture for the ancillary channel £. Using 
this somehow involved notation, we can translate the Gregoratti and Werner 
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theorem stating that a channel £ admits a random-unitary representation (^J) 
if and only if there exists a rank-one POVM {|ai)(ai|a} such that 

£'{\ai){ai\a) = Vila, (7) 

for all 2, and for some fixed probability distribution pi. In fact 

Pi = Tt[£{p) \ai){ai\a] = Tr[p £'{\ai) {ai\a)], Mp^i. (8) 

In other words, the POVM {|aj)(aj|a} is mapped to a classical dice, namely, 
the POVM {pila}- Notice that the cardinality of the POVM {\ai){ai\a] 
coincides with the cardinality in the random- unitary decomposition (Q). 

3 Extremal rank-one POVM's 

Let us now suppose that > (dim=^)^. Then we know that such a POVM 
is non extremal ^Sl El El E] and it can be convexly decomposed into 
extremal components 

\a,){a,\a = \P^ + {l-\)Q^. (9) 

(In the above equation we considered a convex combination of just two 
extremal terms; the general case does not change the conclusions.) Since 
{|Q;j)(Q!j|a} is rank-one and < A < 1, the only possibility to satisfy Eq. Q 
is that the non- null elements of {Pi} and {Qi} are all proportional to the 
corresponding element of {|aj)(ai|a}- Hence, by linearity, also the non-null 
elements of {Pi} and {Qi} are mapped by £' to something proportional to 
la. The normalization is granted by the normalization of the map £' . This 
means that at the end we found two other rank-one POVM's, that is {Pi} and 
{Qi}-, that are both extremal, and hence both with cardinality less or equal 
to {dim J^aY, achieving two other random-unitary Kraus representations for 
the channel £. 

On the other hand, the normalization condition |«i)(ai|a = la rules 
out the possibility that < dim J^a- A von Neumann rank-one measure- 
ment, with {ai\aj) = Sij, achieves the lower bound N = dim J^. 

4 The result 

By now, we showed that a random-unitary channel always admits a random- 
unitary decomposition (Q) involving at most {dimJ^aY unitary operators. 
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We can now tighten this bound by choosing dim M'a as small as possible. The 
smallest^ achievable dim M'a for a given channel E coincides with the number 
of Kraus elements in an orthogonal — or canonical — Kraus representation, 
that is, E{p) = '^jKjpKj with Tr[A'jK;] oc 6ji. Such a number is precisely 
the rank of the Choi- Jamiolkowski ^1 120] positive operator Rs in one-to-one 
correspondence with the channel S and defined as 

R£ = {S®I)\Q){Q\, (10) 

where I is the identity channel, and \Q) is a non normalized = d) 

maximally entangled vector in Jif^Jf. An orthogonal Kraus representation 
of S corresponds then to a diagonalization of Rg. 
Thus, we have the main result 

Theorem A random unitary channel £ always admits a random-unitary 
Kraus representation 

K 

£{p) = Y.V.U.pUl (11) 

1=1 

with 

ranki?^ <K< (ranki?^)^ □ (12) 

The bound ()12|) holds regardless of the dimension d of the input Hilbert 
space It is then reasonable that it fails in accurately describing the 

peculiar case of qubits {d = 2). In fact, it is known that all bistochastic 
qubit channels are actually Pauli channels (see, for example, Ref. 0), that 
is, they can always be written as (apart from an overall rotation of the whole 
Bloch sphere) 

^{p)= X] pmp(yi, (13) 

where {ctq = 1, a^., fXy, a^} are the usual 2x2 Pauli unitary matrices, and Pi 
is a probability distribution. Notice that Tr[o"iCT"j] oc 5ij: This means that the 
Pauli form of qubit bistochastic channels is a diagonalization of the channel 
itself and the equality 

K = ranki?£ (14) 

^See Ref. ^Sl for a detailed analysis of the ancillary space dimension, that is, the 
ancillary resources, needed to implement various possible unitary realizations of a given 
quantum channel. 
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holds in this case. However, already for d = 3 there exist bistochastic chan- 
nels that cannot be diagonalized on unitary operators (for an explicit exam- 
ple, see Ref. |H]). This evidence clearly does not prove our bound to be tight. 
It nonetheless shows that things, already for d = 3, acquire highly non-trivial 
geometric properties and get much more complicated. In all these cases, the 
bound given in the Theorem could be tight. 

As an immediate consequence of the Theorem, it stems the following 

Corollary The minimum amount of classical information needed to be 
extracted from the environment in order to perfectly correct a random-unitary 
channel £ is upper hounded as 



Moreover, since rankRs < d^, the following quite loose — yet independent of 
the particular channel — bound holds 



5 Concluding remark 

It is noteworthy that we need no more than (ranki?^:)^ rank-one POVM 
elements in order to extract all the "useful" classical information from the 
ancilla. This is analogous to what happens in the case of optimal accessible 
information extraction: as proved by Davies [21] , one never needs more that 
d^ rank-one POVM elements in order to extract the maximum achievable 
accessible information from a d- dimensional system. In the case of accessible 
information extraction, Davies' bound seems to be tight, in the sense that 
examples can be constructed in which the maximum information gathering 
is achieved only by a POVM with maximum number of elements [22] ■ If the 
analogy is correct, the bound in the Theorem could be proved to be tight as 
well, while we expect that the bound given in the Corollary can be refined. 
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H{pi) < 21og(rank/2^). □ 



(15) 



H{pi) <4\ogd. 
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